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Qh' Abstract. By means of numerical methods, we solve the equations of motion for the 
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collapse of a shell of baryonic matter, made of galaxies and substructure falling into the 
central regions of a cluster of galaxies, taking into account the effect of the dynamical 
friction. The parameters on which the dynamical friction mainly depends are: the peaks' 



height, the number of peaks inside a protocluster multiplied by the correlation function 
' evaluated at the origin, the filtering radius and the nucleus radius of the protocluster of 
galaxies. We show how the collapse time r of the shell depends on these parameters. We 
give a formula that links the dynamical friction coefficient t] to the parameters mentioned 
above and an analytic relation between the collapse time and rj. Finally, we obtain an an- 
alytical relation between r and the mean overdensity S within the shell. All the analytical 
relations that we find are in excellent agreement with the numerical integration. 
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1. Introduction 

The problem of the formation and evolution of clusters of galaxies has been one of the 
crucial topics of the last years (see e.g. Ryden & Gunn 1987, Colafrancesco et al. 1989, 
Antonuccio-Delogu 1992, Kaiser 1993, Colafrancesco & Vittorio 1993, Croft & Efstathiou 
1994, Dutta & Spergel 1994, Mosconi et al. 1994, Sutherland & Dalton 1994, Efstathiou 
1994 and Colafrancesco et al. 1995). It is well known that the formation of cosmic struc- 
tures is strictly related to the evolution of the density perturbations: in the present 
paradigm of structure formation, it is generally assumed that cosmic structures of size 
~ R form preferentially around the local maxima of the primordial density field, once it 
is smoothed on the filtering scale Rf. These linear density fluctuations eventually evolve 
towards the nonlinear regime under the action of gravitational instability; they detach 
from the Hubble flow at turn around epoch tm, given by: 



3tt 



■(1 + 5) 



1/2 



{1 + zf/' (1) 



32Gpb 

where is the mean background density, z is the redshift and S is the mean over- 
density within the nonlinear region. After the turn around epoch, the fluctuations start 
to recoUapse when their overdensity defined by 

.(X) . (2) 

Pb 

reaches the value 6=1. 
The evolution of the density fluctuations is described by the power spectrum, given by: 

V{k) ^ {\S^\'} (3) 

with: 

5k = y" d'^ke-'^'^Six) (4) 

Since the density fleld depends on the power spectrum, which in turn depends on the 
matter that dominates the universe, the mean characteristics of the cosmic structures 
depend on the assumed model. In this context the most successful model is the biased 

Cold Dark Matter (hereafter CDM) (see e.g. Kolb & Turner 1990; Peebles 1993; Liddle 
& Lyth 1993) based on a scale invariant spectrum of density fluctuations growing under 
gravitational instability. In such scenario the formation of the structures occurs through 
a " bottom up " mechanism. A simple model that describes the collapse of adensity per- 
turbation is that by Gunn & Gott (1972, hereafter GG72). The main assumptions of this 
model are: (a) the symmetry of the collapse is spherical; (b) the matter distribution is 
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uniform in that region of space where the density exceeds the background density; (c) 
no tidal interaction exists with the external density perturbations and (d) there is no 
substructure (collapsed objects having sizes less than that of main perturbation). 

Point (d) is in contradiction to the predictions of CDM models. It is well known that 
in a CDM Universe, an abundant production of substructures during the evolution of 
the fluctuations is predicted. 

The problem of the substructures in a CDM Universe and his consequences on struc- 
ture formation have been widely studied in previous papers (see e.g. Antonuccio-Delogu 
1992, hereafter A92; Antonuccio-Delogu & Atrio-Barandela 1992, hereafter AA92; 
Antonuccio-Delogu & Colafrancesco 1994, hereafter AC94, Del Popolo et al. 1996, Gam- 
bera 1997 and Del Popolo & Gambera 1997 hereafter DG97). Being a rather recent topic 
in cosmology much work has still to be done. 

The presence of substructure is very important for the dynamics of collapsing shells of 
baryonic matter made of galaxies and substructure of lO^M© lO^M©, falling into the 
central regions of a cluster of galaxies. As shown by Chandrasekhar & von Neumann 
(1942, hereafter CvN42; 1943), in the presence of substructure it is necessary to modify 
the equation of motion: 

df^ r^t) ^ ' 

since the graininess of mass distribution in the system induces dynamical friction that 

introduces a frictional force term. 

Adopting the notation of GG72 (sec also their Eqs. 6 and 8) and remembering that 
Tco/2 is the collapse time of a shell of baryonic matter in the absence of dynamical 
friction (GG72), one can write: 

- 1/2 

-fco — 777^ TTTT y^) 

Hi{pi - PciY'^ 

where pd is the critical density at a time U and pi is the average density inside Ti at 
U. The equation of motion of a shell of baryonic matter in presence of dynamical friction 
(Kandrup 1980, hereafter K80; Kashlinsky 1986 and AC94), using the dimensionless 
time variable r = t^- , can be written in the form: 

d^a A-KGpci{l+Ji) 2 rj. da 
d^ - ^) " "^^'^d^ 

where 6i is the overdensity within r,, r] is the coefficient of dynamical friction and 

a{ri,t) is the expansion parameter of the shell (see GG72 Eq. 6), that can be written as: 

a{n,t) = '-^ (8) 
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Supposing that there is no correlation among random force and their derivatives, we 
have: 

^ Jd^FWiF)F^T{F) 
^ ~ 2(«2) 

(K80), where T{F) is the average "duration" of a random force impulse of magnitude 
F, W{F) is the probability distribution of stochastic force (which for a clustered system 
is given in Eq. 37 of AA92). 

DG97 solved Eq. (|) numerically, showing qualitatively how the expansion parameter 
a(T) depends on the dynamical friction coefficient and how r changes in the presence of 
dynamical friction, but without undertaking a more complete study of the dependence 
on the parameters. 

The plan of this paper is as follows. In Sect. 2 we show how r depends on the peaks' 
height i^c, on the parameter S that we define there as the correlation function at the origin 
multiplied by the total number of peaks inside a protocluster, on the filtering radius Rf 
and on the nucleus radius of the protocluster tq. A more detailed description of these 
parameters is given below. In Sect. 3 we give an analytical relation between the dynamical 
friction coefficient and the collapse time: 

T = fliv) (10) 

In Sect. 4 we give a semi-analytical relation that links the dynamical friction coefficient 
with the parameters on which it depends: 

V ^ f2{'^c,Rf,rQ,E) (11) 

Linking these two relations we also find the dependence of the dimensionless collapse 
time T on the parameters used: 

T = /io/2 = /i[/2(i/e,i?/,ro,S)] (12) 
In Sect. 5 we give a semi-analytical relation between r and 6 and rj: 

r^h(S,v) (13) 

Finally, in Sect. 6 we summarize our results and comment on their possible implica- 
tions. 

2. The collapse time 

DG97 showed how the expansion parameter a(T) depends on the dynamical friction, 
solving Eq. (|^) by means of a numerical method but not taking into account the param- 
eters on which r] depends. 
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Here we examine how the dynamical friction coefficient 77 varies according to the param- 
eters and how the collapse time depends on them. We consider Eq. (||) and the functions 
T{F) and W{F), the former given by Chandrasekhar (1943), the latter (for clustered 
system) given by the so-called Holtsmark law (Holtsmark 1919): 



2F 

W{F) = — / dkksm{kF)exp 
= 47r|F| VF(|F|) 



3/2' 



(14) 



where W{F)dF is the probability for a test particle of experiencing a force in the range 
F^F + dF, N is the total number of particles, rrityp is a typical particle mass and rtyp is 
a typical distance among the particles. 
According to K80: 

1 



W{\F\)^ 



with: 



27r2|F| 



dfcfcsin(fc|F|)^/(fc) 



Af{k) ^ lim Ajv(k) = A„or exp(jr) 



(15) 



(16) 



The original expression for given by K80 has been modified by AA92 to take into 
account clustering, and turns out to be given by: 



Auncl{k) + Acl{k) 



2A' 



2,cl 



(17) 



where Anor, Aci{k), Auncii^), A'^ ^iid Sci(^) ^^'S given in AA92 respectively 

by the Eqs. (21), (29), (31), (32) and (36). 

Here, we want to remind that Sc;(fc) is a linear function of the correlation function ^(r) 
and that the general expression for Yjci{k), adopting the notation of AA92 is given by: 



Scz(k) 



d Yi I d Y2 exp I i 



■ exp I i^Uh:^!! ) rc;(ri)rc((r2)C(|ri - rzl) 



(18) 



where masses are measured in units of solar mass Mq and distances in h^^ Mpc, so 
that k will be measured in unit of /i^^(Mpc)2/GM0, with niav the average mass of the 
substructure. Since we have mav> 10^ Mq >> IMq for all the cases that we consider in 
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this paper, we can adopt the asymptotic expansion of Sci demonstrated in the appendix 
of AA92: 

^ci{k)^m (19) 

(see AA92 for details). Ecj does not depend on ^(r), but only on ^(0) (that is ^(r) 
calculated at the origin). 

We solve Eq. (|^) and the other equations related for a outskirts shell of baryonic matter 
with 6 = 0.01 inside the spherical regions (protoclustcr), for different values of Vc, Rf, rg 
and S (where S = Ntot ■ f (0)), the latter quantity being better defined in Sect. 4. 
After having determined rj solving numerically Eq. (^, we get r as a function of i^c, Rf, 
rg and S solving Eq. (^. We perform these calculations for different set of values of i^c, 
Rf, rg and S inside the following intervals: 



1.2 <Vc< 3.2 
0.5 /i"^Mpc < ro < 10 /i"^Mpc 
10"^ /i-^Mpc <Rf <l /i^^Mpc 
20 < S < 2 • 10^ 



(20) 



The results that we have obtained are shown in Figs. 1h-4. Before commenting upon the 
figures, we want to remark that the dependence of r on (5 is qualitatively shown in Fig. 5 
by AC94. We observe that for 5 > 10^^ the collapse time in the presence of a dynamical 
friction is always larger than in the unperturbed case but the magnitude of the deviation 
is negligible for larger 6, whilst for 5 < 10^^ the deviations increase steeply with lower 
S. Then, having considered 6 — 0.01, the estimation we get for r in Sect. 4 must be 
considered as a lower limit. 

In Fig. 1 we show the collapse time in the presence of dynamical friction, versus the peaks' 
height, for different values of S. In this picture, we show how r grows for larger values 
of h'c and for larger values of S. Similarly, in Fig. 2 we note how t increases for larger 
values of Vc and of Rf. The slope of the curves confirm our prevision on the behaviour 
of the collapse of a shell of baryonic matter falling into the central regions of a cluster 
of galaxies in the presence of dynamical friction: the dynamical friction slows down the 
collapse (as DG97 had already shown) and the effect, as we are showing in Figs. 1 and 
2, increases as S, Rf, Vc grow. 
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Fig. 1. Collapse time t of a shell of matter made of galaxies and substructure when dynamical 
friction is taken into account, versus Vc We assume a nucleus radius of ro = l/i~^Mpc and a 
filtering radius i?/ = 0.74/i~^Mpc. Open circles: H = 10'^; filled circles: H = 10^. 

Here we want to remind that we are considering only the peaks of the local density field 
with central height v larger than a critical threshold Vc- This latter quantity is chosen 
to be the threshold at which Vp^ak (y > Vc) << lav where rpeak is the typical size of the 
peaks and lav is the average peak separation (see also Bardeen et al. 1986). 
In Figs. 3 and 4 we show how the collapse time varies with the nucleus radius of the 
protocluster tq. Note how r grows as ro decreases: the smaller the nucleus of the proto- 
cluster, the larger the time of collapse in the presence of dynamical friction; besides we 
show how this effect increases for larger values of both i^c and S. 



3. How the collapse time depends on the dynamical friction 

Like DG97, we have solved Eq. numerically and calculated the collapse time r for 
different values of the dynamical friction cocfficent rj. However, DG97 showed the depen- 
dence of r on 77 qualitatively, whilst we have found also an analytic relation r = fi^r]) 
that links t with ?; in the range of values < rj < 3.1 (see also Gambera 1997). The 
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Fig. 2. Collapse time t of a shell of matter made of galaxies and substructure when dynamical 
friction is taken into account, versus Vc We assume a nucleus radius of ro = l/i~^Mpc and a 
fixed correlation H = 10^. Open circles: Rf = Mpc; filled circles: Rf = 0.65h^^ Mpc; 

crosses: Rf = 0.55h~^ Mpc. 

results of our calculations are shown in Fig. 5, where we report the collapse time r versus 
the dynamical friction coefficient. We have determined the function r = through a 

least squares fit method. Our best solution is given by: 



(21) 



with: 



ai = 0.00595923, 
02 = 0.0493460, 
aa = 0.492164, 

Eq. and our results for < 77 < 3.1 from the numerical integration of Eq. are 
in excellent agreement. As a matter of fact, the correlation coefficient value obtained is: 
= 1.000. 
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Fig. 3. Collapse time t versus ro. We assume a filtering radius Rf — 0.74/i ^ Mpc. and a total 
number of peaks of substructure H = 10^. Open circles: Vc = 3; crosses: Vc = 2; filled circles: 
Va = 1.24. 

4. How the collapse time depends on the parameters 

As already shown by AC94, it is possible to write down the dynamical friction coefRcient 
T] as the sum of two terms: 



V = Vo + Vci 



(22) 



where rjo is the coefficient of dynamical friction of an unclustered distribution of field 
particles whilst r^c; takes into account the effect of clustering. We rewrite Eq. ( p2|) as: 



V = Vo + Vol = Vo [ ^ -\ 



(23) 



where, as demonstrated by A92, the ratio ^ depends only on ^ and vq whilst r]Q is given 
by A92: 



3_33 [GH p{^,. )c,f' Xog[lM{n)li'RW) 



,3/2 



(24) 
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Fig. 4. Collapse time t of a shell of matter made of galaxies and substructure when dynamical 
friction is taken into account, versus ro We assume a filtering radius Rf — 0.74h~^ Mpc and a 
peaks' height uc = 3. Open circles: H = 10^; filled circles: H = 10'^. 



and, for a fixed value of S, depends only on Vc and Rf. In the previous formula, Cp is 
defined in Bower (1991), whilst p{j, Vc), {m)av and {n)av are defined in AC94. Therefore, 



we rewrite Eq. (23) as 



V = Voii'c, Rf) [1 + ^(ro, 5) 



(25) 



where the parameters on which rj depends are the following: 

— i/c is the peaks' height; 

— Rf is the filtering radius; 

— tq is the parameter of the power-law density profile. A theoretical work (Ryden 
1988) suggests that the density profile inside a protogalactic dark matter halo, before 
relaxation and baryonic infall, can be approximated by a power-law: 



p{r) 



P 



(26) 



where p « 1.6 on a protogalactic scale. 
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Fig. 5. Collapse time r of a shell of matter made of gala^xies and substructure when dynamical 
friction is taken into account, versus rj. 

— S is the product Ntof^{0) where Ntot is the total number of peaks inside a protoclus- 
ter and ^(0) is the correlation function calculated in r = 0. AA92 have demonstrated 

that in the hypothesis niav 3> lAf©, where niav is the average mass of the subpeaks, 
the dependence of the dynamical friction coefficient on N^ot and ^ (r) may be expressed 
as a dependence on a single parameter that we define as: 



s = Ntot ■ m 



(27) 



The analytic relation r] = /2(fc, ^i/, t"o, S), that links the dynamical friction coefBcient 
with the parameters on which it depends, can be rewritten as the product of two functions: 



V = M^c, Rf,ro, E) = f!^{Pc, Rf) ■ fHiro, S) 



(28) 



where 



(29) 



and 
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1 I ' Cl ell ( "\ 



(30) 



With a least square method, we find the best function r\ = /2(t'c, i?/, J'oi ■^)- First, we 
find an analytical relation between the dynamical friction coefficient in the absence of 
clustering 770 and the parameters Rf and v. We obtain: 

770 = 61 + + 63^?/ + hi^vRf + h^v^ 



bi = 1.00655066, 62 = -0.01224778, 
bs = 0.03836632, 64 = -0.03337465, 
65 = 0.00288929, be = -0.02622193, 
67 = 0.04374888, 63 = 0.00690435, 
69 = 0.00860054, 

We perform a test between the values of 770 obtained from the last equation and the 
values of 770 calculated integrating Eq. (p^). Our result for the range 1.2 < Vc < 3.2 is 
excellent: ~ 10"*. 

We do the same job for the quantity /2'('"o, S), finding: 

/2 = Cl + C2ro + C3S + CiEro + c^rl + ceSrQ (32) 



(31) 



with: 



with: 



Cl = 2.26516888, C2 -0.09679056 



C3 = 0.00054051, C4 = -0.0001255044, 



C5 = 0.00850948, ce = 0.000009336723 



An analogue test gives ~ 10 * for the range 20 < ^ < 2 • 10"^. 

The function 7/ = /2(i^, i?/, ro, S) is given by the product of Eqs. ( pl| ) and (p2[). These 



contain 54 terms. 



77 = f2{i',Rf,ro,E) 



= (5i + 62/^ + bsRf + b^vRf + 651/^ + bg-R/ + ^yi^^/ 
+ b^v^R) + bgiy^Rj) ■ (ci + C2ro + C3S 
+ c^Ero +C5rQ+ ceEr^) 



(33) 
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However, we have also found an empirical formula with only 13 terms that represents 
a good approximation. We performed a test between the results obtained from the 
13-term equation and the results obtained from the numerical integration. The result is 
= 2.17 • 10^^. The same test performed on the 54-term equation gives « 10~®. Our 
13-term equation reads as: 

rj « di + (121^ + d^Rf + d4ro + d^E + d^vRf 

(34) 

+ diii^Rj + dui^'^Rf + disr^E 



with: 

di = 2.280000, d2 = -0.027743, 
da = 0.086906, d^ = -0.097425, 
ds = 5.440520, de = -0.075599, 
dr = 0.006545, 4 = -0.059397, 
dg = 0.008565, diQ = 1.26 • 10-^ 
dii = 0.099099, di2 = 0.015640, 
di3 = 9.39 • 10-6 

The function t — fi o f2 is given by: 

/i o /2 = oi [/2(j^, i?/,ro,5)]^ + 02 [,f2{'^,Rf,ro,E)] + 03 (35) 
that is: 



.0 fi+^Y 


1 


^0 (^n 




+ 02 




. V mJ. 









where the values of a„ are given in Sect. 2 and the function /2 is given by Eq. ( |33| ) or 
by Eq. (H). 



5. A semi-analytical relation between r and 6 

Our aim is to find a semi-analytical relation r = /a ((5, 77) for the intervals lO^"* < 5 < 10"^ 
and < 77 < 3.1. 

The first step is the determination of a function r = ^((5) for a fixed value of 77. We solve 
Eq. (0) for lO^"' < <5 < 10^^ and r\ = 0.01 and by a least square method we find the 
function t = g{6): 

T = eiS^ + 62^ + 63 (37) 
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with: 

ei = 5678.65, 

62 = 45.3429, 

63 = -0.01591, 

The value of the correlation coefficient between this function and the numerical inte- 
gration is = 1.000. With this method we find r = g{S) for several values of rj inside 
the interval < < 3.1. We can write the dimensionless collapse time as the product: 

T = g(S) ■ h{r^) (38) 

where the function h{ri) can be written as 

h{7j) =K-[1 + /(ry)] (39) 

and, for lO"'^ < ^ < 10"^ ^nd < t] < 3.1, K is constant: K = -1494.7705898. So: 

T = F(S,rj) = K-[l + f{^)]-g(S) (40) 

A test performed between the values obtained for t from the numerical integration 
and the values obtained from Eq. (^ gives the result: ^ 10~^. 



6. Conclusions and discussion 

In Sect. 2 of this work we showed in a quantitative way how the collapse time r of a 
shell of baryonic matter made of galaxies and substructure depends on some parameters. 
When one of the parameters S or i?/ or increases, the collapse time grows. It means 
that the effects of the presence of dynamical friction should be more evident in the outer 
regions of rich clusters of galaxies. Besides, we show how the collapse time of an infalling 
shell increases with decreasing values of tq, and becomes very large for tq < 2h^^ Mpc 
(see Fig. 4). As a consequence, the slowing down of the collapse of an outer shell within 
a cluster of galaxies owing to the dynamical friction is more remarkable in the clusters 
with nucleus of little dimension. 

In Sect. 3 of this paper we give an analytic relation that links the dimensionless collapse 
time T with the coefficient of dynamical friction rj. This relation is in excellent agreement 
with the numerical integration of Eq. (^ for < 77 < 3.1 (see Fig. 5). Then, we find an an- 
alytic relation between 77 and the parameters on which it depends for the value 6 = 10^^ 
(we remind that the effects of the dynamical friction are negligible for S > 10^^). This 
is Eq. (|2^), that can be considered as a "low order" approximation to a more realistic 
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situation of an outer shell of a cluster of galaxies with S < 10^^. We also find an empirical 
formula, Eq. (Q), that is a good approximation of Eq. (|3^. Moreover, the dependence 
of the dimensionless collapse time on Vc, fQ, S and Rj is shown (Eq. (|35|)). 
Finally, we give an analytical relation that links r with t] and S. Here, we wish to stress 
the usefulness of an analytic relation like t — K ■ [1 + {fi o /2)] • g{S). This is a powerful 
tool to estimate the effect of the dynamical friction in the outer regions of clusters of 
galaxies (Gambera et al. in preparation) and to compare the observational data with the 
theoretical ones. This is a good method to test how important is the role of the dynamical 
friction in the collapse of the clusters of galaxies. 
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